Wigner-Dyson statistics for a class of integrable models 
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We construct an ensemble of second-quantized Hamiltonians with two bosonic degrees of freedom, 
whose members display with probability one GOE or GUE statistics. Nevertheless, these Hamil- 
tonians have a second integral of motion, namely the boson number, and thus are integrable. To 
construct this ensemble we use some "reverse engineering" starting from the fact that n-bosons in a 
two-level system with random interactions have an integrable classical limit by the old Heisenberg 
association of boson operators to actions and angles. By choosing an n-body random interaction 
and degenerate levels we end up with GOE or GUE Hamiltonians. Ergodicity of these ensembles 
completes the example. 

PACS numbers: 05.45.Mt, 05.30.Jp, 31.15.Gy 



Recently, there has been considerable interest in 
spinless n-boson systems with fc-body random in- 
teractions 01 0, both for degenerate and non- 
degenerate 0, Q single-particle levels. In particular, 
anomalous statistics for the two-level system have been 
understood from the fact that these systems are inte- 
grable J6|. The two-level ensemble corresponds to a 
time-independent two degrees of freedom system, in the 
sense that creation and annihilation operators for the two 
single-particle levels are canonical operators for the sys- 
tem. The second integral of motion corresponds to the 
conservation of the number of bosons. The Hamiltonian 
is a function of the creation and annihilation operators, 
whose precise form depends on the question whether we 
have two-body interactions or more complicated many- 
body interactions. The number of (random) coefficients 
in this model increases quadratically with the rank k of 
the interaction. A classical Hamiltonian can formally be 
written for any number of particles but in fact only for 
large particle number the quantum problem reaches the 
classical limit 0, • Thus the boson number plays the 
role of an action. For fixed particle number, the model 
can be reduced to a Hamiltonian of one degree of freedom 
with the number of bosons appearing as parameter. 

In this paper we shall consider a different large-n limit 
for the two-level system. Our model consists in choosing 
the rank of the interaction equal to the particle number, 
k — n. Hence the Hamiltonian changes for different val- 
ues of the particle number. This choice, by definition 
of the ensemble, directly leads to a GOE or GUE ma- 
trix . We then consider the limit of large particle num- 
ber. Clearly, we will have an ensemble of systems which 
have GOE or GUE statistics with probability one [To| . 
while being integrable in a well-defined limit of large ac- 
tions. This represents an important caveat concerning 
the idea that Wigner-Dyson statistics are characteristic 
of classical chaotic systems. Note that this has no bear- 
ing on the fairly-well established quantum chaos con- 
jecture [Hill El Q E E E3, which establishes that 
classical chaos implies typically Wigner-Dyson statistics. 



In a similar vein, Wu et al. [Tsl ] have performed the fol- 
lowing calculation: after having chosen an unfolded fixed 
spectrum of given length N generated via the diagonal- 
ization of a matrix taken from the GOE, they fit a one- 
dimensional potential, the spectrum of which coincides 
with the above spectrum for its first N values. The main 
difference between their work and ours is, that we display 
explicitly an ensemble of integrable Hamiltonians having 
a GOE spectrum. 

The Hamiltonian in second quantization is 
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Here, creates a spinless boson in single-particle level j 
(j = 1,2), bj destroys it, and the fc-body matrix elements 
v s j corresponds to a GOE or a GUE 3]. The number 
operator n — b\b\ 4- b\t>2 commutes with the Hamilto- 
nian independently of the rank k of the interaction. The 
combinatorial factors A/" s ,t = [s\{k — s)\t\(k — t)!] 1 / 2 in 
Eq. (TJ3) are introduced in order to have exactly a GOE or 
GUE Hamiltonian for n = k, where n is the number of 
bosons P, ■ For the two-level system, the dimension of 
Hilbert space is N = n + 1. 

The classical Hamiltonian is obtained as follows |a.ll9|: 
H is symmetrized with respect to the ordering of the cre- 
ation and annihilation operators. This permits a correct 
assignment of the zero-point energy and therefore a one- 
to-one classical-quantum energy comparison. Then, we 
use Heisenberg semiclassical rules Q 
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Formally, the classical Hamiltonian can be written as 
H = Hq + V. The unperturbed Hamiltonian Hq, which 
includes a constant associated with the zero-point energy 
correction, depends only upon the actions I\ and 1%; the 
residual interaction V depends also on the angles 4>i and 
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These terms are explicitly given by 
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F = ^Af cos[(s-t)(0i -0 2 )] 



x ( Pt [l 1 -^, s } + p t [I 1 -l,t}) 

x (P fe _ s [/ 2 - I fc - s ] + P fe _ s [/ 2 - I - t]).(4) 

In these equations (details of their derivation will be 
given elsewhere), Pt[I, s] is a polynomial of degree t on I 
defined as (s > t > 0) 



Pt[I,s] = l[[I-(s 
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FIG. 1: Phase-space structure of a specific realization of the 
reduced Hamiltonian for k = n = 30. 



Here we have assumed real matrix elements v Sj t, which 
are independent random variables, Gaussian distributed 
with zero mean and variance given by v Stt v s ^t' = $ss> 8tt> + 
$st' $ts' i with the overline denoting ensemble average. 
Hence the /c-body interaction matrices form a GOE. The 
analogous case for the GUE is obtained by replacing in 
Eq. Q v a< t by |u s ,t| and introducing a phase w Sj t in the 
cosine function. 

The Hamiltonian is written as a polynomial of the 
two actions with random coefficients, and cosines of the 
difference of the two angles. The invariant is given by 
K = Ii + 1- 2 = n+l which reflects the translational sym- 
metry of the interaction with respect to the angles. Again 
we can test the corresponding Poisson brackets 6] . This 
form displays explicitly that we deal with an integrable 
two-degrees of freedom problem. If we fix the invariant, 
i.e. the particle number n, we can reduce the Hamilto- 
nian to one degree of freedom. The quantum Hamilto- 
nian acts on a Hilbert space of dimension n + I, which 
is the number of ways we can distribute the n bosons in 
the two single-particle levels. Formally, this corresponds 
to the so-called "polyads" for algebraic Hamiltonians in 
molecular systems |Tj|. The reduced Hamiltonian is ob- 
tained by performing the canonical transformation de- 
fined by the generating function W = K(f>i + J(4>2 — 4>i)- 
The new actions K and J and their corresponding canon- 
ically conjugated angles \ an d V' are related to the old 
variables by 

h =K-J, I 2 = J, (6a) 
X = 0i, ip = h - <j>i- (6b) 

Substituting the new variables in Eqs. Q and (@J, and 
fixing K yields explicitly the reduced Hamiltonian. 

We now have established the main point of the pa- 
per: if we consider k = n, by construction the quantum 
Hamiltonian coincides with a GOE 0, 0]. Each mem- 
ber of the ensemble is associated, in the limit of large 



n(= k), to a classical Hamiltonian with one effective de- 
gree of freedom. This defines the ensemble of integrable 
Hamiltonians. 

Let us now qualitatively understand why this hap- 
pens. For this purpose we investigate the structure of the 
phase space of a typical member of the ensemble for fixed 
ti>l. A Poincare surface of the two degrees of freedom 
system constructed for constant K — n + 1 is equivalent 
to a contour plot of the reduced one-dimensional Hamil- 
tonian. Fixing K makes for a compact phase space in 
the reduced coordinates ip and J defined above. The 
volume corresponds to n + 1 as we normalized to one 
quantum state per unit cell. The number of random co- 
efficients entering in the Hamiltonian grows quadratically 
on k. The equations of motion then contain polynomi- 
als in J and K — J (of degree up to k = n) multiplied 
by sines or cosines of multiples of ip. In this case, the 
degree of the polynomials sets ~ n 2 as the upper limit 
for the total number of fixed points (stable and unsta- 
ble). While a more precise estimate for their number is 
difficult to establish, our numerical results show that the 
fixed points proliferate more rapidly than n, and come 
close to the upper bound. The typical phase-space por- 
trait will therefore consist of elliptic islands surrounded 
by separatrices in a complicated mesh. In Fig. ^ we il- 
lustrate the phase-space structure of a member of the 
ensemble for k = n = 30. The periodic orbits whose 
action satisfy the EBK quantization condition 
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define implicitly the z-th energy level Ei. Here, ni is an 
integer and a>i is the Maslov index of the orbit. Note 
that there is a priori no monotonic dependence of the 
value of the action as a function of energy for the reduced 
Hamiltonian. That is, for a given energy two or more 
distinct invariant tori may exist in different regions of 
phase space; each of these tori may have a different value 
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of the corresponding action. 

As mentioned above, the number of stable and un- 
stable fixed points proliferates very rapidly for growing 
k = n. This is a consequence of the many-body charac- 
ter of the interaction V in Eq. J2J. In particular, in the 
limit we consider, this number grows much faster than 
n ~ h , this implies that the phase-space volume sur- 
rounding each stable fixed point shrinks for n = k — ► oo. 

For one degree of freedom systems in the large n limit 
with k fixed the spectrum is constructed from sequences 
of levels obtained by torus quantization around the ellip- 
tic fixed points, i.e. most periodic orbits which satisfy 
the EBK condition can be uniquely assigned to one ellip- 
tic fixed point. The spectrum is thus a superposition of 
picket-fence (equidistant) sequences of levels as reported 
in Q for small values of k. For fixed (sufficiently large) 
values of k we expect a Poisson spectrum to arise on short 
length scales. This expectation is based on the fact that 
WKB theory yields a superposition of many picket fence 
(equidistant) spectra with different spacings. 

In contrast, in the present case (k = n) we have ~ n 2 
elliptic fixed points and only n + 1 energy levels. There- 
fore, a stable periodic orbit that fulfills the EBK con- 
dition (JJJ, instead of surrounding an elliptic point, has 
to accommodate and explore more extended regions in 
phase space, coming close to the separatrix associated 
with many different unstable fixed points. Stated in a dif- 
ferent way, while the periodic orbits of the Hamiltonian 
are strictly stable, wave packets started on initial condi- 
tions separated by a distance of the order h 1 ^ 2 ~ n" 1 / 2 
will sample very different regions of phase space. This 
coarse-graining will effectively mimic unstable motion. 
This fact makes plausible that we do not obtain a smooth 
spectrum that can be unfolded to yield (superpositions 
of) picket-fence spectra. Indeed, we find Random Ma- 



trix spectral fluctuations. Note that the results on the 
non-ergodic behavior in the dense limit (k fixed) of the 
bosonic fc-body embedded ensembles Q cease to apply 
in this case by the same line of reasoning, and therefore 
we recover the standard ergodicity of the GOE |icj . 

Finally, we consider the extension to the GUE case. 
We recall that in action and angle variables of the har- 
monic oscillator, time-reversal invariance is tested by the 
invariance of the Hamiltonian under the transformation 
(j)j — » —<pj (for j = 1,2). As mentioned above, when 
considering complex fc-body matrix elements, Eq. Q re- 
quires certain modifications; namely, v St t is replaced by 
its modulus \v Syt \ and the phases uj s t are introduced in 
the argument of the cosine functions. The introduction 
of the phases u> St t in the cosine functions of Eq. Q breaks 
the time-reversal invariance. 

Summarizing, we have defined an ensemble of inte- 
grable Hamiltonians of two degrees of freedom, or equiv- 
alently, an ensemble of one degree of freedom Hamilto- 
nians by fixing the number of bosons n. In the limit 
n = k — > oo, each member of the (reduced) ensemble 
is mapped exactly by the quantization onto a member 
of the GOE or GUE. Therefore, the fluctuation proper- 
ties of its eigenvalues follow the predictions of Random 
Matrix Theory. This result can be interpreted in two 
ways: First, the limiting Hamiltonian is considered to be 
a classical one, or second, such a limit is not accepted as 
a classical one. In either case, for this family of Hamil- 
tonians Wigner-Dyson fluctuations does not imply chaos 
in classical dynamics. 

We acknowledge financial support by the projects 
DGAPA (UNAM) IN-109000, IN-112200 and CONA- 
CyT 32173-E. LB was supported by the PSPA-DGAPA 
(UNAM) program and thanks O. Bohigas for discussions 
and the kind hospitality at the LPTMS (Orsay). 



[1] V.K.B. Kota, Phys. Rep. 347, 223 (2001). 

[2] L. Benet and H.A. Weidenmuller, J. Phys. A: Math. Gen. 

36, 3569 (2003). 
[3] T. Asaga et al., Eur. Phys. Lett. 56, 340 (2001); Ann. 

Phys. (N.Y.) 298, 229 (2002). 
[4] K. Patel, M.S. Desai, V. Potbhare and V.K.B. Kota, 

Phys. Lett. A 275, 329 (2000) 
[5] N.D. Chavda, V. Potbhare and V.K.B. Kota, Phys. Lett. 

A 311, 331 (2003). 
[6] L. Benet, C. Jung and F. Leyvraz, J. Phys. A: Math. 

Gen. 36, L217 (2003). 
[7] W. Heisenberg, Z. Phys. 33, 879 (1925). 
[8] L.G. Yaffe, Rev. Mod. Phys. 54, 407 (1982). 
[9] ML. Mehta, Random matrices (Academic Press, New 

York, 1991); T. Guhr, A. Muller-Groehling, and H.A. 

Weidenmuller, Phys. Rep. 299, 189 (1998). 
[10] J.B. French, P.A. Mello and A. Pandey, Phys. Lett. B 

80, 17 (1978). 

[11] O. Bohigas, M.-J. Giannoni and C. Schmit, Phys. Rev. 
Lett. 52, 1 (1984). 



[12] S.W. MacDonald and A.N. Kaufmann, Phys. Rev. Lett. 
42, 1189 (1979). 

[13] G. Casati, F. Valz-Gris and I. Guarneri, Lett. Nouvo 
Cimento 28, 279 (1980). 

[14] M.V. Berry, Ann. Phys. (N.Y.) 131, 163 (1981); "Struc- 
tures in semiclassical spectra: a question of scale" in The 
Wave-Particle Dualism, eds. S Diner et al, (D. Reidel, 
Dordrecht, 1984), pp 231. 

[15] M.V. Berry, Proc. Roy. Soc. London A 400, 229 (1985). 

[16] A.V. Andreev et al., Phys. Rev. Lett. 76, 3947 (1995). 

[17] F. Leyvraz and T.H. Seligman, in Proceedings of the 
Fourth Wigner Symposium, eds. N. Atakishiev, T. H. 
Seligman, and K. B. Wolf (World Scientific, Singapore, 
1996), pp. 429. 

[18] H. Wu, M. Valieres, D.H. Feng and D.W.L. Sprung, Phys. 

Rev. A42, 1027 (1990) 
[19] MP. Jacobson, C. Jung, H.S. Taylor, and RW. Field, J. 

Chem. Phys. Ill, 600 (1999). 



